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1. Introduction 

Among the many beautiful connections between statistical physics and combinatorics, 
those involving integrable lattice models are probably the deepest ones. A major advance 
on the enumeration of various symmetry classes of alternating sign matrices (ASM) was 
accomplished by Kuperberg |l[] by relating these notoriously difficult combinatorial prob- 
lems to partition functions of the six-vertex model on a square grid with various boundary 
conditions. Alternatively, Okada reinterpreted these partition functions as characters 
of various classical Lie groups. 

In this note, we address a related combinatorial mystery, the Razumov-Stroganov 
conjecture |J, that identifies the suitably normalized groundstate vector entries of the 
0(1) loop model on a semi-infinite cylinder of perimeter 2n to numbers of configurations 
of the fully-packed loop model on an n x n square grid. A weaker version of this conjecture, 
the Razumov-Stroganov sum rule H identifies the sum of these entries with the number 
A n of ASMs of size n x n, and has been proved in in an inhomogeneous version 
(incorporating position-dependent spectral parameters Zi) making extensive use of the 
integrability of the loop model. It allows for the identification between the sum of entries 
of the fully inhomogeneous 0(1) loop model and the very same partition function as 
that used by Kuperberg for the enumeration of ASMs, the so-called Izergin-Korepin (IK) 
determinant ||. This proof was adapted to the case of the inhomogeneous 0(1) crossing 
loop model as well |7J], and allowed to identify the entries of the groundstate vector with 
the multidegrees of some matrix varieties ||. The case of open boundaries for the loop 
models was also treated in || both for non-crossing and crossing loop models respectively, 
giving rise to other sum rules. In particular, this established the open Razumov-Stroganov 
sum rule, identifying the sum of suitably normalized groundstate vector entries of the 0(1) 
loop model on a semi-infinite strip of width 2n with the number Ay(2n + 1) of vertically 
symmetric alternating sign matrices (VSASM) of size (2n + 1) x (2n + 1). 

In the context of (/-deformations of quantum hall effect wave functions |[L0| , Pasquier 
recovered the Razumov-Stroganov sum rule (for q = — e^) by finding the minimal polyno- 
mial solution of the quantum Knizhnik-Zamolodchikov (qKZ) equation |yj for U q (sl(2)) 



at level 1. In [12], solving qKZ equations for higher rank algebras U q (si(k)) at level 1 
led to a nice generalization of the Razumov-Stroganov sum rule for q = — e k + 1 , iden- 
tifying a weighted sum of entries of the groundstate vector of the inhomogeneous Ak-i 
interaction-round-a-face (IRF) model on a semi-infinite cylinder of perimeter kn with par- 
ticular Schur functions generalizing the IK determinant and Okada's result. In particular, 
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in the homogeneous limit where all Z{ — ► 1, the sum rule was found to be proportional to 
the number: 

^ir !n nf +1>i+ ' )! (i.i) 

reducing to the number of ASMs for k = 2. 



In the present paper, we extend the results of |12| to the case of A^-i IRF models 
with open boundaries. To this purpose, we introduce in Sect. 2 a boundary version of the 
level 1 qKZ equation which incorporates a parameter r, playing the role of a boundary 
magnetic field, and such that the ordinary qKZ equation is recovered when r — > 0. This 
extends the boundary qKZ equation of |13| and |L4]], defined respectively for U q (sl(2)) and 



U q (sl(k > 2)), both at level 0. We then use the path representation of |T^] for the quotients 
of the Hecke algebra associated to the A^-i models to find minimal polynomial solutions 
to these boundary qKZ equations in Sect. 3. At the generalized Razumov-Stroganov point 
q = — e k + 1 , we compute a weighted sum of entries of this polynomial solution, identified 
with the groundstate entries of the corresponding inhomogeneous A^-i IRF model on a 
semi-infinite strip of width kn, with open boundaries (Sect.4.1). This sum rule produces a 



simple generalization of the Schur functions of |T2[ , in a form very similar to the character 
formulas of Okada for symplectic groups. In the homogeneous limit where all z± — > 1, all 
these polynomials collapse to polynomials of r with integer coefficients, and lead at r = 1 
to a natural generalization of the number of VSASMs: 

{ } ~ l\ «* + mm + i)u + n) - 1)! {LI) 

which reduces to Ay(2n + 1) for k = 2. 

We also investigate the other "combinatorial point" q = — 1, corresponding to the 
rational limit of the underlying integrable model (Sect. 4. 2). There, in the homogeneous 
limit, we find that all entries of our solution tend to integers, and discuss their possible 
interpretation as degrees of nilpotent matrix varieties. 

2. Quantum Knizhnik-Zamolodchikov equations without and with a boundary 

2.1. R matrix and path representation 

Throughout this note, we shall use the standard abstract trigonometric solution of 
the Yang-Baxter equation which reads 

^ / x q~ 1 z — qw T z — w 

Ri, i+1 (z,w =^tt (2.1) 

q L w — qz q L w — qz 



where the e^, i = 1, 2, iV — 1 are the generators of the Hecke algebra Hn(t), namely 
subject to the relations 



= re* 



(2.2) 



with the parametrization 



r = — 



(9 + 9 ') 



(2.3) 



A standard pictorial representation for R reads 




i+1 



i+1 



(2.4) 



z 



where we have also indicated the action on the spaces i and i + 

More precisely, we let R act on the so-called path representation of the A^-i IRF 
models. To this purpose, the are further constrained to generate the U q (sl(k)) (also 
called Ak-i or SU(k)) quotient H^\t) of Hn(t) obtained by imposing the extra van- 
ishing conditions Yfc(ej, ej + i, . . . , ei + k-i) = 0, i = 1,2,...,N— k, where the Young 
g-symmetrizers are defined inductively by Yi(ei) = ej and F m+ i(ej, . . . , ej +m ) = 
Y m (ei, . . . , ei +m _i)(e i+m - fJ, m )Y m (ei, . . . , e i+m _i), where /x m = C/ m _i(r)/C/ m (r), £7 the 
Chebyshev polynomials subject to t/ m (2cosx) = sin(m + l)x/ sinx. Note that k = 2 cor- 
responds to the Temperley-Lieb algebra TL^v(r). The path representation is only briefly 
sketched here and will be described in detail in |15| . The representation basis is indexed by 
closed paths of length iV = nk, from and to the origin, on the oriented-link Weyl chamber 
of SU(k), generated by the k vectors u\ = Ui, ui = 1V2 — u>i, Uk-i = ojk-i — oJk-2, and 
Uk = — Wfc-i, in terms of the fundamental weights uj{. The paths are allowed to visit only 
points A = ^i^i with all Aj > 0, and may be indexed by their sequence of steps, and 
by a slight abuse of notation we write it = (iiiz . . . in) for the sequence , Ui 2 , . . . , Ui N ; 
we also note 7r m = i m , the m-th step of the path n. For instance, the path of length 
nk closest to the origin is (12.../c) n , namely n repetitions of the sequence 1, 2, . . . , k, 
and we denote it by ixj. Likewise, the path farthest from the origin is (l) n (2) n . . . (k) n , 
and we denote it by ttq. A useful notation consists in representing each step j by a unit 
segment forming an angle of ^(fc+i)"^ with the horizontal direction: each 71 becomes 
a broken line touching the x axis at its ends and staying above it. There are exactly 

c nk = (kn)\ rio<j<fc-i J 'V( n + j)- such paths. The number is the dimension of the 

— — (2) 
path representation and generalizes the Catalan numbers c n = c 2j [ ■ The next step consists 
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in decomposing each path in the broken line representation into "tiles" made of (possi- 
bly glued) lozenges Ljj with edges corresponding to pairs l<i<j<kof steps. To 
each tile we then associate a specific polynomial of the e's. For instance, each single 
lozenge with steps at positions i,i + 1 along the paths corresponds to e^. The tile de- 
composition of each path is then unique up to elementary moves between tiles allowed by 

(k) 

the algebra relations, and we may therefore associate a unique element of H N to each 
path 7r, thus forming a representation basis |tt). Actually the map is onto the left ideal 

Hfi'Cl, O = Y" fc _i(ei, ...,e fc _i)Y fc _ 1 (e fc+1 , e 2 k-x)—Yk-l{^(n-i)k+ii e nfc -i), isomorphic 

( k) 

to H\r • For illustration, let us detail the case k = 3, n = 2. There are five paths, respec- 
tively decomposed as follows: defining i^j+i = Y 2 (e{, e, + i) = e^e^ie^ — e{, Z{j = e^e.,- — 1, 
and O = 11,2^4,5, we have 



M = (X) = 




and the ej read in this basis: 
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(2.6) 



Let us now simply list a few basic properties of this representation, of crucial impor- 
tance for the following: 
(PI) ei \tt) = t \tt) if 7Ti < 7Ti_|_i (n locally convex) 
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(P2) Ei |tt) = Cj,7r,7r' I 71 "') if > ^i+i locally flat or concave), for some Ci^^i G {0, 1} 
(P3) If Ci^y = 1 then tt' is locally convex between steps i and i + 1, namely 7r^ < 7r^ +1 
(P4) If Ci^^i = 1 then either (71^, 7^+1) = (7^+1, ttQ and Tv m = TT' m for all m 7^ z, z + 1, i.e. 
7r' exceeds tt by the unit lozenge L Wi+1)iri in the broken line representation, or tt' C tt, 
namely the broken line representation of tt' lies below that of tt (or equivalently tt' 
may be completed into tt upon adding a number of lozenges). 

We conclude with another important property of this representation. There is a duality 
between the SU(k) and SU(n) representations for iV = kn, confirmed by the identity of 
dimensions = cj^. Indeed, there exists an involution 99 from A^-i to A n _i paths 
of same length kn, defined as follows. First of all, ip takes ttq of Ak-i to tt f of A n -\. 
The bijection is then constructed iteratively by simultaneously subtracting (resp. adding) 
lozenges at the same positions from the A^-i (resp. to the A n _i) paths, so that tt is 
locally convex iff </ ? ( 7r ) is locally concave. A simpler way of picturing this involution uses 
yet another expression of Ak-i paths, as rectangular standard Young tableaux with k row 
of n boxes, in which successive steps along the path are recorded by their positions in the 
corresponding row of the tableau, i.e. if the m-th step is j, then we write the number m 
in the first available box from the left in the j-th row. In this language, ip is simply the 
transposition of the tableaux. Moreover, this duality extends to the path representations, 
namely 

r {n) _ Mk) , 7 \ 

In practice, this relation is extremely useful: it allows for instance to express the represen- 
tation ( p.5| )-( [2!6|) of H§ (r) entirely in terms of the standard Dyck path representation of 
TL 6 (r). 



Once we pick this path representation, the i?-matrix (|2.1|) may be interpreted as 
encoding face weights of some IRF model based on SU(k), with edge degrees of freedom 
in the set of steps, and acting locally at steps i, i + 1 of the paths. This is the natural 
generalization of the 0(n) loop model in the Dyck path representation, corresponding to 
k = 2 and n = r. 
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Fig. 1: The regular- (a) and boundary- (b) qKZ equations. Each oriented 
line carries a spectral parameter. Each intersection between two such lines 
corresponds to a R operator. On the second line, we have represented the 
diagonal operators of respectively multiplication of the spectral parameter 
by s (shift), inversion of the spectral parameter with weig ht r (K( r \z)) and 
weight rs (K~( rs \z)), repectively pictured as a crossed dot, empty dot, filled 
dot. The former is nothing but the composition of the two latter. 

2.2. qKZ equations 



(2.8) 



The qKZ equation [£EI| reads 

, Zj-i, Zj, Zj + i, 

Q = R j+ljj (zj +1 , szj) . ..R N j(z N , szj)Sj(z)aRij(zi, Zj) . 

with the clear pictorial representation of Fig.l-(a). Here the operator a rotates cyclically 
the spectral parameters and the "shift" operator Sj(z) acts as the identity on the space j 
while it multiplies the spectral parameter z by some scalar s (it is represented by a crossed 
dot in Fig.l-(a)). In the pictorial representation of Fig.l-(a), each intersection between two 
oriented lines stands for an R operator ( |2.4| ). Note that the "return" Zj line (from left to 
right) passes below the other (vertical) lines, hence there are no intersections, i.e. no R ma- 
trices involved: the only effect is the operator a, which rotates the spaces and their spectral 



parameters according to (zj, z\, 



■1 z j- 



1; z j+l: 



ZN 



Zl, 



-1, Zj + l, 



Z N , Zj). 



In the presence of boundaries, the qKZ equation is modified as depicted in Fig.l-(b). 
The "return line" now crosses the vertical ones, and the shift of spectral parameter is 
now replaced by two boundary operators at the far left and right of the return line. The 
corresponding qKZ equation now reads 



Zj-i, SZj,Z j+1 , ...,Zn) = Q^( z l, ...,Zj-i,Zj, Zj + i, ...,z N ) 

Q = Rj+ij(zj+i, szj) . . . Rnj(zn, szj)K- ( — )-Rj,jv( — 5 z n) ■ ■ ■ Rj,j+x{ — > z j+i 



Zj Zj Zj 



(2.9) 



r , ^ , r 



x R i,i-i(—i z i-i) ■ ■■Rj,i{^^ z i) K i { z j) R i,j{ z i> z j) ■ .-Rj-ijizj-uZj 



z j z j 
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The boundary operator (z) acts as the identity on the space j, and changes the spectral 
parameter z into r/z. Together with R, it satisfies the boundary Yang-Baxter equation 
of Ref. |T(J. In Fig. lb, we have represented respectively and K^ rs ^ by an empty and 
a filled dot. Here r is an additional free parameter, which allows to recover in the limit 
r — > the qKZ equation without boundary (|2.8|) from that with boundaries (|2.9|) . Indeed, 
in this limit, the spectral parameter along the return line tends to 0: consequently, all 
R matrices along the return line degenerate into pure under crossings, and we recover the 
picture of Fig.l-(a), upon noting that S(z) = K^ rs \r / z)K^ r \z) is nothing but the shift 
operator that multiplies the spectral parameter z by s. 

The equation ( |2.8|) is equivalently standardly reduced to a system of equations for \P, 
namely 

Ti^(zi, ...,z N ) = Ri ti+ i(zi + i,Zi)ifr(zi, ...,z N ), i = 1, 2, . . . ,N - 1 

(2.10) 

aV(z 2 , z N , szi) = c*(zi, z N ) 

where T{ acts on functions of the z's by interchanging Z{ and c is an irrelevant constant, 
and a now acts as a rotation of states in the representation picked for the R's. Similarly, 
we may reduce (|2]^) to the following system: 

Ti\P(zi, ...,z N ) = • • • , zjv), i = 1, 2, . . . , N - 1 (2.11a) 

r 

— , z 2 , zjv-i, zjv) = ci(zi)^(zi, z N ) (2.116) 

^1 

#(zi, z 2 , zjv-i, — ) = c N (z N )V(zi, ^tv) (2.11c) 

Z N 

for some functions of a single spectral parameter ci and cat. Note that Eqs.( [2.116, c| ) 
express the covariance of ^ under the action of K~[ r \zi) and K^ s \zn)- The latter are 



obtained by showing that the operator Q in ( |2.9|) is intertwined by K± and K^ s \ as a 
direct consequence of the boundary Yang-Baxter equation. The boundary qKZ equation 
( [2 has the advantage of being ridden of the cyclicity condition (second line of (|2.10|) ), 



replaced by reflective boundary conditions, somewhat easier to deal with. We now restrict 
ourselves to vectors \& in the path representation above, namely look for \I/ = ty n \tt), 
7t running over the set of Ak-i paths of length N = kn. When written in components, 
Eqs. (|2.11] ) become 

-qz i )dj$ ir (zi,...,z N ) = ^ Cj^'^rf (zi,---, z N ), 1 < % < N - {2.12a) 
r 

*tt(— j«2,—»«JV-i,*Ar) = c 1 (z 1 )^ 7T (z 1 , ...,z N ) (2.126) 

Zl 

T S 

^(zi, z 2 , zn-u — ) = c N (z N )^ n (z 1 , z N ) (2.12c) 
zn 

where the divided difference operator c\ acts on functions of the z's as 

dif = Tif ~ f (2.13) 

Zi — Zi+i 



We also denote by ti = (q 1 Zi+i — qzi)d%- The equations (|2.12a|) express nothing but 
= (e, -r)*. 

By analogy with the ordinary qKZ case, when dealing with the Ak-i path represen- 
tation above, we write 

8 = q 2{h+l) (2.14) 



and call / the level of the (boundary) qKZ equation ( |2.11| ). The parameter r is free, and 
has the dimension of the square of a spectral parameter. In the next section, we construct 
polynomial solutions of the system ( |2.11j ) for level 1 = 1. 



3. Polynomial solutions for the level 1 boundary qKZ equation 

We now look for polynomial solutions of level 1 boundary qKZ equations ( |2.11|) with 
minimal degree. 
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3.1. Solutions 

We have the remarkable property for all k, that the equations ( |2.12a|) may all be 
solved for in terms solely of ^t^. This relies on the property (P4) above of the path 
representation. Indeed, we may order the paths by inclusion. Then each path tx locally 
convex at i,i + 1 may be viewed as the addition of a lozenge L 7I -- 7r . +1 to a smaller path 
7r~, such that 7r~ = 7r m for m 7^ i, i + 1 and (tx~,tx~ +1 ) = (TXi + i,TXi). Using property 
(P4), we see that eiTX~ = tx + ^tx", where the sum runs over paths tx" C tx~ . Turning 
now to Eq.( p.l2a|) , we see that, as all the tx' involved in the sum on the r.h.s. are such 
that Ci j7r ' i7r = 1, only one of them is strictly included in tx, namely tx~ , while all others 
contain tx. Stripping iteratively lozenges from paths, we may therefore express each ty^- 
as t^Tr — X] 7 r"^ f 7r"' with already known tx" containing tx. This allows for expressing 
iteratively every ty n in terms of VP^, in a triangular manner w.r.t. inclusion of paths. 

The "fundamental" component ^J/^ is fixed by (i) some highest weight conditions 
(ti + T)ty no = for all i not multiple of n (expressing that txq is locally flat except at the 
convex points multiple of n, see below), (ii) the reflective boundary conditions (|2.12b, q 1 ) , 
and finally the requirement that it be of smallest possible degree. 



As a direct consequence of Eq. (|2.12a|) , ^ vanishes whenever Zj+i = q 2 Zi and tx is 



concave or flat at i, i + 1, namely TXi > Tx i+ i. Indeed, from property (P3), the vector \tx) 
cannot be the image of anything under e^; but precisely when Zi+i = q 2 Zi, the matrix 
R becomes proportional to (c.f. the definition ( |2.1| )), hence ^ n must vanish. This is 
straightforwardly extended to any portion of tx flat or concave say between positions i 
and j, upon considering the appropriate product or R matrices: if z\ = q 2 Zk for some 
i < k < I < j, then ^/^ vanishes. Applying this to txq which is flat everywhere except 
at points multiple of n where it is convex, we find that ty^o must factor out the terms 
n(m-i)n+i<i<i<mn('? 2 ^ ~ z j) for a11 m = 1,2, k. Alternatively, wherever tx q is locally 
flat, i.e. for all i not multiple of n, we have (ti + t)\I/ 7I - = 0, as |7To) has no antecedent 
under e^: this "highest weight condition" is analogous to that of ||12|| . It implies that once 
the above factors are taken out of the remaining is a polynomial symmetric in the 
variables {z( m -i) n +i, ^(m-i)n+2j — » z mn } separately for all m = l,2,...,k. Finally, by an 
inductive argument similar to that used in Ref. , we can show that factors out terms 
(q 2m r — ZkZi) for all pairs /c, / such that (m — l)n + l < k < I < run, m = 1,2, k. Indeed, 
the left boundary condition (|2.12fe|) allows to show that Y\2<j<n( ( l 2r ~ z i z j) divides ^no- 
Next, we let z\ — > oo, namely pick in ^ the coefficient of top degree in z\\ in this case the 
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space 1 decouples from the picture of Fig. 1(a), as its crossings with the j line are undone 
and replaced by undercrossings, and we are left with the same boundary problem, but with 
line 1 erased, and now a reflective boundary condition with parameter r in terms of Zi- 
Repeating iteratively this process, namely taking successively Z2, £3, z n -\ — > 00 exhausts 
all the announced extra factors for m = 1. Assuming each variable occurs in ^ with 
maximal degree 2(n — 1), we find that the leftover dependence on z n must be proportional 
to and we may write ^(00, 00, z n , z n+1 , z N ) = Cz^~ 1 F(z n+1 , z N ). The 

full vector \l/(oo, 00, z n , z n+ \, zn) now obeys a left reflective boundary condition of 
the form 

r 

*(oo, 00, — , z n+ i, z N ) = ci(z n )^(oo, 00, z n , z n+ i, z N ) (3.1) 

Zn 

We now note that ttq = e n ir~ , for the unique tt~ such that tv~ = 7T; for i 7^ n, n + 1, and 
(7r~, 7r~ +1 ) = (7r n+ i, 7r n ). According to (|2.12a| ), we have t n ^ 7To = ^V-. Expressing that 



the component obeys the new left reflection condition (|3.1|) solely in terms of \E' 7ro 
yields the following relation for F: 

cx(zn) q ~ 1Zn ^ - QZn {{z n+1 ) n - l F{z n , ...) - {z n ) n - x F{z n+1 , ...)) 

Zn Z n -\-l 

q- 1 z n Zn+i - qr . \n-lvf r ' r N n_1 



qr — z n z n +i 



((^n+ir- 1 ^-,.-.)-!-) F( Zn+1 ,...)) 

(3.2) 



2 

Picking z n = q~ 2 z n +i finally yields j-, ...) = (q ,2 r/z^ +1 ) n_1 F(2 n+ i, ...) and therefore 

^ttj, now has all the announced extra factors for m = 2 (everything works now as if the 
left boundary parameter had been rescaled as r ^ rq 2 ). Proceeding iteratively, we finally 
get all the announced factors for m = 1,2, k. We claim that the solution of minimal 
degree has no extra factors in its namely that 

k 

= n n ^ - ^)N 2m - ztzj) (3.3) 

m=l (m — l)n+l<i<j<mn 



and that Eqs.(|2.12"D are satisfied for s = q 2 ( k+1 ^ and 



(>) . ^^=(^2— ( 3 - 4 ) 



This fixes the solution entirely. 
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As in the case without boundaries of [12], the solution may alternatively be fixed by 



some global requirement that it vanishes for some specific choices of the spectral parame- 
ters. It was found indeed that in the case without boundary \I/ must vanish for any ordered 
sequence of k + 1 spectral parameters of the form {z, q 2 z, q 2h z}, the so-called "quantum 
incompressibility" condition, first noticed by Pasquier |10| in the case k = 2, and allowing 



for interpreting ^ as (/-deformed wave function of fractional quantum Hall effect. In the 
present case with boundaries, we find that \1/ must vanish for any ordered k + 1-tuple of 
spectral parameters of the form 

{^V-V^X^,..,^^} (3.5) 

z z z 

for m = 2, 3, k + 1. This generalizes the quantum incompressibility condition, in that 
the wave function must vanish whenever k + 1 electrons come into contact with one another 
or now with their reflections. 

As in the case of |jl2| , we may also derive recursion relations for the components, 
relating size N = kn to size N' = k(n — 1). This is done by simply taking z\ = z,Z2 = 
q 2 z,...,Zk = q 2 ( k ~^z, in which case only paths with 7ii = i for i = 1,2, ...,k have non- 
vanishing components in The rest of the path, tv' = (iTh+i, — > ttjv) is an path of 
length N' = N — k, and we have 

( N \ 

i $i V (z 1 ,...,z N ) = I Yl (q 2h z - Zj)(rq 2 - zzj) \ ^ n ,(z k +i, ...,z N ) 

2l =g-2 22 = ... =(? -2(fc-l) 2fe=2; I j=k+l J 

(3.6) 

where the prefactor is fixed by the degree and vanishing requirements. 
3.2. Other sizes 

As already touched-upon in the previous section, there is an easy way of extending 
our polynomial solutions to all sizes N, not necessarily a multiple of k. Starting from a 
solution at = kn, we may indeed send successively the first j spectral parameters to 
(or equivalently to oo), leaving us with only non- vanishing components of \& for paths 7r 
starting with the (unique) convex sequence 1,2,3, This provides us with a solution 
to the boundary qKZ equation of zn, hence in size N — j. Note that now in the 

corresponding representation the paths start at the point u\ + U2 + ■■■ + Uj in the Weyl 
chamber of SU{k) and end at its origin. This construction is valid for all j = 1, 2, k — 1 
hence allows to exhaust all possible sizes. 
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4. Sum rules 

4-1. Generalizations of the Razumov-Stroganov sum rule at q — —e k + 1 



We proceed as in |fL2|| . We introduce a covector v such that 

e iV = TV, i = l,2,...,N-l (4.1) 

This covector is unique under the condition that all its entries be strictly positive and 
that v nf = 1, and exists only if q = —e k + 1 or its conjugate, namely r = 2cos-^-. The 
entries of v are directly constructed from the path representation and may be expressed 
as products of ratios of the Chebyshev polynomials U m (r) with m < k, associated to the 
tiles in the decomposition of the paths, and all of which are positive for the above choice of 
r. For illustration, for k = 3, n = 2, we take r = \/2, and associate the factor U\(t) = y/2 
to every lozenge and the factor ^(t) = 1 to every tile made of two lozenges in the path 
decomposition (|2.5|), corresponding respectively to the left eigenvalues of v on and on 
Zij, resulting in v = {1, v^2, 1, 1, v^2} in the path basis. 

The condition ([O]) immediately implies that vRi^+x (z, w) = v, readily checked from 
the definition (|2.1|) . Henceforth, the combination v ■ ^ is a symmetric polynomial of the 
z s, as a direct consequence of Eq. ( |2.11a|) . Introducing 

I ( k k J(z ll z 2l ... 1 z kn \r) = (-z/g)^- 1 )/ 2 ^-^^,...,^) (4.2) 

with i = \/ — 1, we find that lj® is a polynomial of degree kn(n—l)/2 of r, whose coefficients 
are all symmetric polynomials of the z's, and the top degree and constant coefficients are 
respectively 

l kn | r fen(n-l)/2 



r( fc )| _ I \ 



(4-3) 

J kn \r° = {ziZ2---ZknT S Y [Zl, Z kn ) 

where sy is the Schur function associated to the Young diagram Y with k rows of n — 1 
boxes, k rows of n — 2 boxes, k rows of 1 box. The second line of ( f4.3| ) is nothing but the 



generalized Razumov-Stroganov sum rule of [12 1, recovered in the limit r — > as explained 



above. The polynomial Zk n \ r ) has only non-negative integer coefficients in r and 

the z's. Note also the covariance under z% — > r/z{ for each i, inherited from the boundary 
reflections ( |2.116, c] ) and the overall symmetry, namely 

r fr\ Ul 
lkn ^i-i' — ' •■•! z kn\f) — [ ~2 J ^fet ^i— 15 z j? ^i+ij •••) z k n \r) (4.4) 
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which allows to relate the coefficients of for different powers of r. In particular, taking 
simultaneously all z% — > r/zi relates the top and constant coefficients, and the first line 

fer 



of (|4.3|) follows from the second. The polynomial ij^ (|4.2[) is entirely determined up to 



global normalization by the following properties: 

(i) it is a symmetric polynomial of the z J s quasi-homogeneous under Z{ — > Xzi for all % 
and r — > A 2 r, with total degree 3kn(n — l)/2 and partial degree 2(n — 1) in each of 
the z's and kn(n — l)/2 in r 

(ii) it vanishes for all the choices of ordered A;-tuples of spectral parameters ( p.5|) , with 
q = -efe+i 

and the normalization is fixed by ( |4.3|) . 

Let us first discuss the Temperley-Lieb case k = 2 and n arbitrary, where drastic 
simplifications occur. At the Razumov-Stroganov point q = — e~?, r = 1, we have v n = 1 

(2) (2) 

for all 7T, hence J^n ^ s simply the sum of components of Vl/. The quantity J^n reduces to the 
result of at r = 1. For generic r, the determinantal formulas of 0] are straightforwardly 
generalized to 

r ( 2 ) , I Uij=i ( z i + z * z J+n + Zj+n) (r 2 + r Zl z J+n + zfz] +n ) 



I^(z 1 ,...,z 2n \r) 



Yll<i<j<n( Z i Z j)( r z i z j)( z i+n z j+n)( r z i+n z j+n) 

- - , (4.5) 



x det 



i<i,j<n y zf + ZiZ j+n + z 2 j+n r 2 + rziz j+n + z 2 z 2 +n 



which is an open version of the Izergin-Korepin determinant , and the Pfaffian formula 
becomes 



(2) 2 u (Zj +z l z J +z J )(r +rz l z J +z l z J ) 
I 2n { Zl ,...,z 2n \r) - [[ (z--z-)(r-z-z-) 

l<i<j<2n {Zl Zj)[ ZlZj) 



X Pfx 



<i<j<2n 



(zi - zj)(r - ZiZj) 



(4.6) 



{z 2 + ZiZj + z 2 ){r 2 + rZiZj + zfz 2 ) 



These are proved for instance by using the symmetry and the recursion relations obeyed 
by v ■ *B as a consequence of Eq.( |3.6|) , and showing that the r.h.s. is the unique symmetric 
polynomial obeying the same constraints. In the homogeneous case where all Zi = 1, 

(2) 

I\n reduces to a reciprocal polynomial of r with non-negative integer coefficients, equal 
respectively to 3 n ( n_1 )/ 2 times the number A n of alternating sign matrices (ASM) of size 
n x n at r = (a result proved in |5[), and to 3 n ( n_1 ) times the number Ay{2n + 1) of 
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vertically symmetric alternating sign matrices (VSASM) of size (2n + l) x (2n + l) at r = 1 
(a result proved in 0), with respective values 

, "tt' (3j + , , „ ... "rr' (gj + 4 )'(2j + 2)! .. 

(k) 

For general k, we have the following expression for the sum rule . 

I^izi, ...,z nk \r) = a Y (zi, ...,z nk \r) (4.8) 

where sy are "symplectic" Schur functions depending on the extra parameter r, and defined 

as 

deti<i ,,•<*„ U kn+1 ~ i+j - (r/ziY^-'+A 

s Y (zi, z nk \r) = z 1 z 2 ...z k n -7—. t (4.9) 

det x <i,j<kn 14 - {r/ziYj 

with Y as in (^4.3|) , and li = n — 1 — [(i — l)/fc], is the number of boxes in the i-th 
row, z = 1,2, ...,kn (we have added fc bottom rows of box for convenience). The sum 
rule ( |4.8| ) is proved as usual by checking that the r.h.s. fulfills all requirements of degree 
and vanishings above, and is therefore fixed by uniqueness. This formula is a nice multi- 
parameter generalization of that written by Okada for VSASMs, in terms of specialized 
characters of the symplectic group. 

In the particular case n = 2 and k arbitrary, It is easy to show that (|4.8|) reduces to 

k 

4 k k \z 1 ,...,z 2k \r)= J £r m s Ym (zi,...,z 2k ) (4.10) 

m=0 

expressed in terms of ordinary Schur functions for the young diagrams Y m with k — m rows 
of 2 boxes and k rows of 1 box. Consequently, in the homogeneous limit where all Z{ — > 1, 
we have 

i /c + 1 /2fc + l\ /2fc + l 



7^(1, 1, l|r) = V r™^^- - ^ M - ^ M (4 } 

2fc v ' ' ' 1 ; 2jfe + l\ m \k-m> 

m=0 

and we get 



A 

the /c-th Catalan number, while at r = 1 we obtain 



? ) = ^ T l2 ) (l,l,-,l|0) = c fc (4.12) 



^( fc )f2l - 1 I (k) (l 1 1111-2 ( 4fc + 1 ) ! (aio) 
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1 


2 


3 


4 


5 


i 


1 


1 


1 


1 


1 


2 


1 


3 


13 


68 


399 


3 


1 


26 


1938 


246675 


43475640 


4 


1 


646 


3251625 


46278146640 


1238084585216726 


5 


1 


45885 


61003011480 


444316978958627636 


9080679253196247653297250 



Table I: The numbers Ay (n) for 1 < n, k < 5. 



(k) 

For arbitrary k, the homogeneous case of all z* = 1 yields again for 1^ a reciprocal 
polynomial of r, with non-negative integer coefficients. The constant coefficient of this 
polynomial is the result of the homogeneous limit of the generalized Razumov-Stroganov 
sum rule computed in [I2|, and reads (k + l) n ( n_1 )/ 2 times an integer A^ generalizing 



the number of ASMs A r 



(2) 

An , and which may be rewritten as 
Ti!nti((fc+l)i + i)! 



= n 

i=o 



Uto (kj + i) 



(4.14) 



The numbers ( f4.14| ) are nothing but the properly normalized dimensions of the GL^ n 
representations Y ||17|| . At r = 1 the homogeneous value of lj® reads (k + l) n ( n_1 ) times 
generalizing the number of (2n + 1) x (2n + 1) VSASMs A v {2n + 1). 



a number A 



(*). 
v 



n. 



We find that 



A 



(k) 
v 



(»)=n 



(2j)\(2nk + 2j - 1)1 
l = \ ((k + l) 3 y.((k + l)(j + n)-l)l 



(4.15) 



These numbers are the properly normalized dimensions of the SP2kn representations Y 
17]. The numbers Ay[n) are easily checked to reduce to Ay{2n + 1) for k = 2 and all n, 
and to ( ]4.13| ) for n = 2, and all /c. The numbers Ay (n) are displayed in Table I for n,k = 



1, 2, 5. For illustration, for fc 



3, n 

I (3) _ 



2 we have /, 



(3), 



1,1, 



20+84r+84r 2 +20r 3 , 

from which we read the values of A 2 °> = 20/4 = 5 and A ( y ] (2) = (20+84+84 + 20)/4 2 = 13. 

As mentioned in Sect. 3. 2, we have access to analogous sum rules for arbitrary sizes 
iV not necessarily a multiple of k by sending spectral parameters to and reducing the 
space of path states accordingly. In size nk — j, the sum rule simply reads (|4.8|) , with say 
Z\ = z% = ... = Zj = and where the first j rows have been deleted from Y. This leads to 
more integer numbers in the homogeneous limit. 
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4-2. More combinatorics at q = — 1 

In ||12|] , a connection was established between the "rational limit" q — > — 1, Zi — > 1 of 
the U q (sl(k)) level 1 minimal polynomial solution of the qKZ equation without boundary 
and so-called extended Joseph polynomials, related to the geometry of the variety of upper 
triangular matrices with vanishing k-th power. As a particular corollary of this relation, 
the entries of \& were found to tend to positive integers in the rational homogeneous limit, 
in which one lets first all Zi — > 1, and then q — > — 1 after dividing out by the appropriate 
power of q + 1. These integers were interpreted as the degrees of the components of the 
abovementioned variety. 

In the present case, we may also consider an analogous rational limit. Of course, at 



r = 0, we may reproduce the rational limit of [|12|. If we want a different result from that 



of [pT2[| , we need to let r — > 1 as well. Substituting 

q = - e - ea/2 1 Zl = e - ew \ r = e~ ep (4.16) 

into our minimal polynomial solution ^> of the U q {s\{k)) level 1 boundary qKZ equation, 
and taking e — > 0, we obtain a two-parameter (a, p) family of polynomials of the ic's. 
At a = p = 0, these are some sorts of reflected versions of the Joseph polynomials ||18|| . 



which await some algebro-geometric interpretation. A remarkable property is that for 
a = 1, p = 0, and in the homogeneous limit where Wi — > for all i, the above polynomials 
reduce to non-negative integers, yet to be interpreted as degrees of components of some 
matrix variety. 

Remarkably, for k = 2 and n arbitrary, we have found that for a suitable normalization 
of ^ these would-be degrees sum to the total degree of the "Brauer scheme" considered in 
||, suggesting a possible connection. More precisely, taking r = 1 and Zi = 1 for all i in 
our solution we have 

lim — L__y ^(1,1, ...,1) = det f^ + Zj + A / 417 ) 

( g 2 _ l)2n(n-l) ^ ' ' > ' <ij<n-l \ 2i J K ' 

where the numbers on the r.h.s. are the degree of the Brauer scheme {M e M2 n ((D), M • 
M = 0} where • stands for a suitable deformation of the ordinary matrix product ||, 
also obtained as sum rules in the 0(1) crossing loop model on a semi-infinite cylinder 0. 
For illustration, in the cases n = 1,2,3 we find the following limits V^n when q — > —1 of 
= l| r = l)/(q 2 - l)2n(n-l). 

4 1} = {1}, 4 1) = {5,2}, 4 1) = {149,52,58,40,8} (4.18) 
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expressed in the standard Dyck path basis of TL 6 (r), ordered by inclusion. The compo- 
nents of (|4.18|) sum respectively to 1,7, 307, the degrees of the Brauer scheme for matrices 
of size 2, 4, 6 respectively. For completeness, the corresponding r = solutions of [TI| yield 
the following limits ip^ when q — > — 1 of *&(zi = l|r = 0)/(q 2 — l) n ( n_1 ): 

4 0) = {1}, ^ 0) = {2,1}, 4 0) = {10,4,4,4,1} (4.19) 

The above observation is particularly intriguing, as we already know that another particular 
point (q = — eV) of the same solution \I/ also displays entries leading to non-negative 
integers for both r = and r = 1, as mentioned in the previous section. These entries 
were actually identified via the full Razumov-Stroganov conjecture || to fully-packed loop 
configurations on a n x n square grid (at r = 0) or to vertically symmetric fully-packed 
loop configurations on a (2n + 1) x (2n + 1) square grid (at r = 1) connected in both cases 
according to the same link pattern as that encoded in the Dyck path it. 

The case k > 2 also seems to lead to non-negative integer entries in the rational 
homogeneous limit, provided we take r — 1 and all Z\ = 1, and then consider 

^ (Q 2 _ DMn-i) E L -> !) ( 4 - 2 °) 

These are presumably related to the degrees of higher order generalizations of the Brauer 
loop scheme [[15]. For illustration, we have for k = 3 and n = 2, in the basis ( |2.5| ): 



^6 1} = {60,28,13, 15,6} (4.21) 

whereas the corresponding r = limit (with (q 2 — i) fcn ( n - 1 ) replaced by (q 2 — i) fcn ( n - 1 )/ 2 
in the limit of ( 4.20|) ) yields in the basis 



^ 0) = {6,3,2,2,1} (4.22) 

As mentioned in Sect. 3. 2, we have access to analogous integer ^'s and sum rules for 
arbitrary sizes N not necessarily a multiple of k, presumably relating to generalized Brauer 
schemes of matrices of corresponding sizes. 
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5. Conclusion 



In this note, we have constructed the minimal polynomial solutions to the level 1 
U q (sl(k)) boundary qKZ equation, and investigated its "combinatorial points". We have 
derived some higher rank generalizations of the Razumov-Stroganov sum rules for the open 
boundary 0(1) spin chain. We have actually found a family of solutions depending on one 
parameter r playing the role of boundary magnetic field. The closed and open boundary 
cases of the Razumov-Stroganov sum rules are recovered for r = and 1 respectively. Sum 
rules for the components of the solution a la Razumov-Stroganov are generically obtained 
at a unique point r = 2 cos for each k, where a covector v with the properties ( }4.1j ) 
and with positive entries can be constructed. 

In view of trying to cook up a full conjecture a la Razumov-Stroganov, we may try to 
consider the vector w with entries w w = (—i/q) kn ^ n ~ 1 ^ 2 v 7T '^ 7r : unfortunately, these entries 
are not real for generic r. For illustration, for k = 3, n = 2, the corresponding entries read 
in the basis ( |2.5|) : 



w f = 4(l+9r + 9r 2 + r 3 ), w 3 = 4(1 + 7r + 7r 2 + r 3 ), w 2 = 4(1 + (2 + i)r + (2 - i)r 2 + r 3 ), 
Wl = 4(1 + (2 - i)r + (2 + i)r 2 + r 3 ), w = 4(1 + r + r 2 + r 3 ) 

(5.1) 

and sum to the abovementioned value of ig (1, 1, 1 |r). In this particular case however, 
the entries of w reduce to non-negative integers in both cases r = 0,1, with w^ ' = 
4{1, 1, 1, 1, 1} and = 8{10, 8, 3, 3, 2}. In general, it is tempting to conjecture that w 
has integer entries, and we have indeed observed this fact for r = and the cases n = 2, 
k = 3, 4, 5, and n = 3, k = 3. For r = 1 however we have found a counterexample for n = 2, 
k = 4, in which the 14 entries of w, although positive and summing to 4 2 A { y\2) = 16x68, 
are not all integers. So if a full Razumov-Stroganov conjecture exists it has to be more 
subtle in the r = 1 case. 

Note that, as we have g 2 ( fc+1 ) = 1 at the generalized Razumov-Stroganov point, the 
left and right boundary conditions in ( |2.11|) coincide. When r = 0, we are precisely in the 
case where \1/ is the groundstate vector of the inhomogeneous A^-i IRF model on a semi- 



infinite cylinder of perimeter iV as pointed out in ||12|| , and governed by the cyclic Hamil- 
tonian H = J2i<i<N( r ~ e i)> e N = crejv-io" -1 an additional generator for the cyclic Hecke 
algebra. For r = 1, ^ is the groundstate vector of the inhomogeneous A^-i IRF model 
on a semi-infinite strip of width N, with reflecting boundary conditions, and governed by 
the Hamiltonian H = Xa<i<iv-i( r — e 0- I n both cases, the entries (— z/g) fcri ^ n_1 )/ 2 \E f 7r 
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and v w are manifestly positive, as these are respectively the entries of the right and left 
Perron-Frobenius eigenvectors of the corresponding Hamiltonians. Alternatively, we may 
view the r = 1 model as describing a stochastically "growing interface" like in [|TI5j for the 
Temperley-Lieb case k = 2, upon interpreting the action of Ri t i+i(zi, Zj) (|2.1|) on a path as 
the identity with probability 1 — tij = (zj — q 2 z%) / (q 2 Zj — Zi) and as e^/r with probability 
Uj, which (i) leaves the path unchanged (if tv is convex at 1), (ii) adds a lozenge to 
it (if 7r is concave at i, i + 1) (iii) or shrinks it to a smaller path n' such that Ci j7I - j7r / = 1 (if 
tx is flat or concave at + 1), the last two operations weighted by 1/r. The combination 
of entries V(tt) = v^^j (v ■ \I/) may then be interpreted as the invariant probability for the 
interface to be equal to the path it. These invariant probabilities read in the case k = 3, 
n = 2: {^ 5 ^;^5^5]^}fo r the paths of (|2~5|). The above counterexample to integrality 
simply shows that these invariant probabilities will not be rational numbers in general. 
It is still possible nevertheless that generalizations of the observables considered in [|19| 
could lead to simple rational expectation values. For instance, from our numerical data, 
we conjecture that the probability of convex transitions iti < iri+i in the above interfaces 
reads 

for the Ak-i growth problem with n = 2. This is easily checked for k = 3, as the paths of 
( |2.5|) have respectively {4, 3, 3, 3, 2} convex transitions among their 5 transitions, leading 
to a probability I{4, 3, 3, 3, 2} ■ {^, ^, ^, ^, ^} = f . 

The sum rules obtained in this paper produce nice integer sequences, which we expect 
to be counting generalizations of ASMs and of VSASMs. A first step in understanding what 
exactly is counted by these numbers could be to find some higher k > 2 analogues of the 
Izergin-Korepin determinant (or its open version (|4.5|)) for our sum rules lu^izi, Zk n \ r )- 



kn 

We may also try to generalize the Pfaffian formula (|4.6|) . The hope would be eventually to 
be able to define an vertex model with a particular type of (domain- wall) boundary 

conditions on a finite domain, whose partition function exactly matches (zi, Zk n \r). 
Then the configurations of such a model would form the desired generalizations of ASMs 
or VSASMs. This might then open the route to generalizations of other symmetry classes 
of such objects, according to the symmetries of the above domain, in the spirit of J[J . 

As a side remark, the number of different objects enumerated by the ASM number A n 
is quite considerable: totally symmetric self- complementary plane partitions (TSSCPP), 
descending plane partitions (DPP), rhombus tilings of particular domains of the triangular 
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lattice, osculating walks on a square grid, ice configurations on a square with domain wall 
boundaries etc, and these give rise to many different formulas for A n , as determinants or 
Pfaffians. We may wonder which ones of these objects and/or formulas will be amenable 
to higher k generalizations, maybe all of them? 

The other combinatorial point q = — 1 (r = 2) seems to be of a different nature. It 
corresponds to the rational limit of the above Ak-i IRF models, namely giving rise to 
a rational solution of the Yang-Baxter equation. Another rational solution of a similar 
nature was considered in and led to sum rules for the 0(1) model of crossing loops. All 
these models seem to be related to the geometry of some specific matrix varieties, allowing 
in particular for interpreting in the homogeneous limit the various integer entries of ^ at 
hand as degrees of the components of these varieties. The identification of these varieties 
for the open boundary cases is under way ]15| . 
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